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Abstract 


Nanoconstruction of metals is a significant challenge for the future manufacturing of plas- 
monic devices. Such a technology requires the development of ultra-fast, high-throughput 
and low cost fabrication schemes. Laser processing can be considered as such and can 
potentially represent an unrivalled tool towards the anticipated arrival of modules based 
in metallic nanostructures, with an extra advantage: the ease of scalability. Specifically, 
laser nanostructuring of either thin metal films or ceramic/metal multilayers and compos- 
ites can result on surface or subsurface plasmonic patterns, respectively, with many poten- 
tial applications. In this chapter, the photo-thermal processes involved in surface and 
subsurface nanostructuring are discussed and processes to develop functional plasmonic 
nanostructures with pre-determined morphology are demonstrated. For the subsurface 
plasmonic conformations, the temperature gradients that are developed spatially across 
the metal/dielectric structure during the laser processing can be utilized. For the surface 
plasmonic nanoassembling, the ability to tune the laser's wavelength to either match the 
absorption spectral profile of the metal or to be resonant with the plasma oscillation fre- 
quency can be utilised, i.e. different optical absorption mechanisms that are size-selective 
can be probed. Both processes can serve as a platform for stimulating further progress 
towards the engineering of large-scale plasmonic devices. 


Keywords: laser annealing, laser induce self-assembly, plasmonics, surface plasmon 
resonance, plasmonic nanoparticles 
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1. Introduction 


Plasmonic nanoparticles have become an increasingly common research area as well as becom- 
ing a key component in many important applications, such as solar energy harvesting [1-6], 
chemical sensing via surface-enhanced spectroscopies [7, 8-14], cancer treatment [15, 16] and 
optical encoding of information [17-20] to name a few. The main reason behind their adapt- 
ability, to these and other prominent applications, is their unique optical properties that allow 
for the manipulation of light below the diffraction limit. Among metals, Au and Ag are consid- 
ered more suitable for plasmonic nanoparticles (PNPs) due to their stability when formed into 
nanoparticles and their strong LSPR absorption bands in the visible region of the spectrum 
[21, 22]. Under irradiation, the free electrons in a PNP are driven by the alternating electric 
field and collectively oscillate at a resonant frequency in a phenomenon known as localized 
surface plasmon resonance (LSPR) [23-27]. The LSPR phenomenon induced in PNPs depends 
strongly on their composition, size, shape, geometrical distribution and on the refractive index 
of the employed dielectric environment [28-30]. To this end, efforts are undertaken in order to 
control these structural parameters depending on the growth conditions [22, 31-34]. 


A critical parameter in delivering practical plasmonic devices is the material preparation 
methods, which should allow the production of nanostructures with tunable plasmonic 
properties. So far, nanomaterials and nanodevice manufacturing have traditionally followed 
two distinct routes: (a) the top-down approach, where a process starts from a uniform mate- 
rial and subsequently finer and finer tools are employed to create smaller structures, like 
lithographic processes [35-37] and/or ion beam nanofabrication [38] and (b) the bottom-up 
approach, where smaller components of atomic or molecular dimensions self-assemble 
together, according to a natural physical principle or an externally applied driving force, to 
give rise to larger and more organized systems, like atomic layer deposition [39], cold weld- 
ing [40], flash thermal annealing [41, 42], pattern transfer [43] and template stripping [44-46]. 
Practically, the top-down approach offers unequalled control and reproducibility down to a 
few nanometres in feature size but at high cost for large area manufacturing, while the bot- 
tom-up route applies for macroscopic scale nanopatterning albeit without the fine feature and 
reproducibility control. 


The use of laser radiation can be viewed as an intermediate route for nanotechnology manu- 
facturing. Light-matter interactions offer precise energy delivery and control over the physi- 
cochemical processes in the nanoworld. Amongst the methodologies followed for creating 
nanoparticle arrangements and/or nanocomposite thin films, laser annealing (LA) has been 
proven to be simple and versatile, providing freedom of design, fast processing, compatibility 
with large-scale manufacturing and allowing for the use of inexpensive flexible substrates 
[17, 19, 47], due to the ultra-fast and macroscopically cold process. For example, LA of a thin 
metal film results into nanoparticle arrangements on top of a substrate [11, 13, 14, 47-54], 
while embedded nanoparticles in a dielectric matrix can be fabricated in large scale through 
the LA of a stack of alternating ultrathin layers of a metal and a dielectric [17, 18]. 


One of the key design issues in LA is the spatial distribution of the energy delivery (from the 
laser pulse) and the precise temporal temperature profile emerged in the structure (caused by 
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this energy). Understanding the dependence of the temperature profile on the nanoparticle/ 
substrate geometries can provide valuable insights in order to be able to design the appropri- 
ate LA process and get custom-made plasmonic devices. In the systems of interest, several 
effects such as material absorption, plasmon resonances and interference contribute to how 
the absorption is distributed, while the material thermal properties and the laser pulse profile 
determine the temperature transients developed. Here, the investigations required to extract 
the underline physical understanding of all these photothermal processes, involved in surface 
and sub-surface laser nanostructuring, are summarised. 


2. UV laser annealing experimental set-up 


A typical laser processing apparatus, such as the one used for the experimental part of the 
work that will be presented later, is shown in Figure 1. It comprises three main components: 
(i) the two laser sources, (ii) the laser beam delivery system and (iii) the sample manipula- 
tion stage. 
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Figure 1. (a) Digital photo of the laser annealing system at Nottingham Trent University, (b) schematic representation of 
the beam delivery path and (c) snapshots of the laser spot at different positions of the beam path. 
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(i) The two laser sources are an excimer laser (LAMBDA PHYSIK LPX 305i) which is capable 
of providing unpolarized light pulses of 25 ns in duration and up to 1200 mJ/pulse at 248nm 
(KrF) and an excimer laser (LAMBDA PHYSIK LPX 205i) capable of providing unpolarized 
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light pulses of 20 ns up to 700 mJ/pulse at 193 nm (ArF). The raw beam, in both cases, is 
approximately 30 x 15 mm? with a Gaussian-like profile on the long axis and an almost top-hat 
profile on the short axis. 


(ii) The beam delivery system is composed of a variable attenuator, a beam homogeniser and 
a mask-projection system. The variable attenuator employs a partially reflective first plate 
(angled at various positions) reflecting part of the laser beam and a second plate that com- 
pensates for the parallel displacement of the beam caused by the first plate. The attenuator 
has a variable throughput from 10% up to 90% of the original beam. The beam homogenizer 
(Exitech Ltd., type EX-HS-700D) employs two lens arrays and a condenser lens. The con- 
denser lens provides a laser spot of pre-determined size (13 x 13 mm?) at a certain distance 
with top-hat profiles on both axes (better than 2% uniformity). This homogenized beam coin- 
cides with a mask which is projected onto the sample via a combination of a field and projec- 
tion lenses. 


(iii) Finally, an X-Y-Z stage is in place in order to manipulate the sample and laser process dif- 
ferent areas on it. The translational stage is operated via a controller with specialised software 
also being capable of externally triggering the laser. Thus, the whole operation is becoming 
fully automated. An appropriate mask and a x5 objective lens set the laser spot delivered onto 
the samples to be a 2.5 x 2.5 mm’ square spot. 


The system has been designed in a way to be versatile and flexible, allowing for an extended 
degree of freedom of operation, thus providing the user with a wide spectrum of experimen- 
tal parameters to vary from. These parameters are shown in Figure 2. 
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Figure 2. The nine experimental parameters of laser annealing system. 


1. Fluence (mJ/cm’): The fluence of the laser, i.e. the energy density delivered to the sample, 
is adjusted by the variable attenuator. The maximum fluence may also vary depending on 
the size of the mask and the objective lens used to project the laser spot onto the sample's 
surface and of course from the excimer laser source: KrF laser is capable of providing a 
raw beam of 1200 mJ/pulse, whereas the ArF is capable of providing a raw beam of 700 
mJ/pulse. Regarding the spot size, this is defined by the appropriate masks, on the mask 
station and the combination of field and projection lenses. 
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2. Wavelength (nm): The KrF laser emits light at 248 nm, while the ArF emits light deeper 
in the UV spectral range; at 193 nm. At this deep UV wavelength, most of the energy is 
absorbed by O, molecules in the air, which are then turned into ozone. To avoid the en- 
ergy losses and thus the minimisation of the energy throughout the process in addition to 
avoiding health hazards due to the toxicity of ozone, the beam delivery path of the ArF 
laser is enclosed in sealed tubes, to allow for N, flow during the laser process. 


3. Pulse duration (ns): Both excimer lasers are ultra-short pulsed UV lasers. The KrF has a 
pulse width of 25 ns, while the ArF has a pulse width of 20 ns. 


4. Number of pulses: The applied number of pulses to the sample. 


5. Repetition rate (Hz): The repetition rate is externally controlled and may vary from 1 to 50 
Hz for both laser sources. 


6. Polarization: While an excimer laser is usually producing completely randomly polarized 
light an optical filter that passes light of a specific polarization (polariser) may be used 
to block waves of other polarizations and convert a beam of light of undefined or mixed 
polarization into a beam with well-defined polarization state. 


7. Environment temperature: The home-built system has the capability to have the sam- 
ple sitting at -196°C (liquid nitrogen) or heated at elevated temperatures during the LA 
process. 


8. Environment composition: The pressure cell, where the sample sits (shown in Figure 1), 
is supplied by a series of gases through appropriate pipelines connected to the cell. These 
gases are Ar, Ar:H,, Ar:H,S, N, and O,. Depending on the application, the LA experiment 
may need to be conducted in an inert atmosphere (e.g., avoiding any oxidation occurring 
during the LA in air) or in a reacting atmosphere, where doping of the processed sample 
is seen necessary. 


9. Pressure: In certain cases, low vacuum or high pressure application during LA (of any gas- 
eous environment) is considered advantageous to avoid detrimental effects of ablation [17]. 


3. Sub-surface nanostructuring 


3.1. Pulsed laser deposited nanocomposite AIN:Ag films 


A robust nanoconstruction consisting of a narrow size distribution of PNPs embedded into a 
ceramic dielectric matrix emerges as the ideal configuration for applications such as authen- 
tication patterns and optical recording media specified for extreme conditions. To this end, 
AIN is an excellent candidate as a host matrix for PNPs due to its optical transparency, 
chemical inactivity, refractory character and high hardness. Ag is a characteristic example 
of noble metal that exhibits the strongest LSPR when formed into nanoparticles. However, 
silver's high diffusivity on one hand and the high melting point of AIN on the other hand 
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makes the fabrication of AIN:Ag(PNPs) nanocomposites an unfeasible task: Ag is either out- 
diffused to the surface (upon thermal treatment) or is atomically dissolved into AIN dur- 
ing growth [55]. Such obstacles prevented the extended use of AIN in plasmonics. Herein, 
strategies for the successful growth of AIN:Ag(PNPs) films by LA-based micropatterning 
are discussed. 


3.1.1. Growth and processing of materials 


AIN:Ag nanocomposite films (70-100 nm thick) were grown on n-type c-Si(100) wafers via 
pulsed laser deposition (PLD) in an in-house-built high-vacuum system (P, < 5 x 10° mbar) 
using the third harmonic of a Nd:YAG laser (6 ns pulse duration, fluence 30 J/cm? at 10 Hz 
repetition rate). More details on sample preparation can be found elsewhere [19, 20]. The 
resulting nanocomposite films consisted of an amorphous AIN matrix with embedded pure Ag 
nanospheres, as displayed in the transmission electron microscopy (TEM) image in the inset of 
Figure 3. The embedded nanospheres exhibited LSPR, spectrally located at around 420 nm [20]. 
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Figure 3. Optical reflectivity of AIN:Ag: as grown film. Inset depicts a cross section TEM image of the sample. 
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3.1.2. Laser processing 
3.1.2.1. The effect of number of pulses 


AIN:Ag films were processed with 1, 2 or 5 pulses at 200 mJ/cm? and 400 mJ/cm?, i.e. just above 
the annealing threshold and just below the ablation threshold at 193 and 248 nm laser wave- 
lengths. The pulse repetition rate in any case was 1 Hz, eliminating any possibility for transient 
thermal phenomena and thus allowing for each successive pulse to interact with a perfectly 
relaxed AIN:Ag film at room temperature. The optical reflectivity spectra are presented in 
Figure 4. In the case of the 193 nm wavelength at low fluence (200 mJ/cm?) (Figure 4a), each 
successive pulse leads to a cumulative LSPR enhancement, in terms of the integral strength of 
the LSPR band, indicating that Ag diffusion results in shape and size improvement of the plas- 
monic nanoparticles; the Si substrate's E1 critical point at around 350 nm [56] is also gradually 
emerging with the number of pulses due to the film's transparency, because of the crystalliza- 
tion of AIN as we will discuss in more detail in the following paragraphs. On the contrary, at 
high fluence (400 mJ/cm7’) (Figure 4b), the evolution of the reflectivity spectra with the number 
of laser pulses exhibits severe changes. The effect of one pulse at 400 mJ/cm? to the plasmonic 
response is more pronounced compared to that of five pulses at 200 mJ/cm/?, indicating that 
the structural modifications taking place during the 193 nm LA process are better associated 
with the laser fluence than with the total laser energy delivered. This will be further examined 
in the next section where the effect of laser fluence is discussed. For two applied pulses, the 
LSPR band is substantially broader indicating an exceptionally broad size distribution of Ag 
nanoparticles that can only be explained by means of Ag outdiffusion. Finally, for five pulses 
the reflectivity spectrum is characteristic of a homogeneous Ag layer on top of the AIN:Ag film 
indicating that none of the remaining Ag is in the nanoparticle form. 


3.1.2.2. The effect of laser wavelength 


A major factor affecting the LA process is the laser wavelength, due to the different absorption 
coefficient of the amorphous AIN matrix and Ag. In the case of LA of AIN:Ag, using the 248 nm 
wavelength and fluence of 400 mJ/cm/’, the picture is quite more complex. At the first laser pulse 
(Figure 4c, red solid line), the pre-existing LSPR band is eliminated, possibly by dissolving of 
overheated, unstable Ag PNPs, whose constituent Ag atoms do not have the required diffusion 
time to reform into bigger Ag particles. Bigger particles, however, are observed after the second 
pulse resulting in intense LSPR peak (Figure 4c, green solid line). The lower photon energy at 
248 nm and the fact that the pulse duration is shorter than the time required for Ostwald ripen- 
ing in a cold environment [57, 58] (the process is indeed relatively cold because the 248 nm beam 
is not absorbed by the AIN matrix due to its high fundamental gap) allows for the Ostwald 
ripening process to take place in two steps: dissolution of particles for the first pulse and reform- 
ing of bigger particles close to the surface for the second pulse. Finally, the optical reflectivity 
spectrum for five pulses (Figure 4c, blue solid line) exhibits a broad LSPR band, resembling the 
broadening of the LSPR band for two applied pulses at 193 nm (Figure 4b, green solid line), but 
blue-shifted compared to the band for two pulses of 193 nm indicating a surrounding medium 
that is less optically dense than AIN (e.g. outdiffused Ag nanoparticles in air) [20]. 
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Figure 4. Optical reflectivity spectra for 1, 2 and 5 successive pulses for laser annealing using: (a) 193 nm wavelength, 200 
mJ/cm? fluence, (b) 193 nm wavelength, 400 mJ/cm? fluence and (c) 248 nm wavelength, 400 mJ/cm? fluence. In all cases 
the black solid line, represents the as deposited film's optical response, which is shown for comparison. Cross-section 
TEM images of AIN:Ag/Si: (d) as-grown, (e) laser annealed with one pulse of 193 nm/400 mJ/cm? and (f) laser annealed 
with two pulses of 248 nm/400 mJ/cm’. 


Cross-sectional TEM images of the as-grown AIN:Ag and laser-treated areas (Figure 4d-f) reveal 
significant structural modifications for both wavelengths. The untreated films (Figure 4d) give 
no trace of crystalline AIN [59]. LA at 193 nm (Figure 4e) leads to agglomeration of Ag nanopar- 
ticles and formation of larger nanoparticles close to the surface of the film. This enlargement 
of Ag nanoparticles can be described as an Ostwald ripening phenomenon. Another important 
aspect of the 193 nm LA process is the partial crystallization of AIN into the wurtzite structure 
(w-AIN) around the enlarged nanoparticles (Figure 4e). Similar results occur after the 248 nm 
LA treatment (Figure 4f), although to a lesser extent due to the lower delivered photon energy, 
which is below the AIN's fundamental gap (6.2 eV). As a result, the laser beam interacts solely 
with the Ag nanoparticles. Consequently, the delivered energy to the substrate is higher than 
in the case of 193 nm LA, and thus the AIN/Si interface is significantly damaged. By the same 
token, the crystallization of the matrix is less likely to occur due to the lower effective tempera- 
ture rise during the LA process. 


To get an insight into the light absorption dynamics inside the film during LA, finite-differ- 
ence time-domain (FDTD) simulations of a 60 nm thick nanostructured AIN:Ag film were 
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performed, where Ag consists of 3 nm diameter inclusions at 15% volume filling ratio ran- 
domly distributed inside the AIN matrix [20]. Laser light, for both wavelengths (193 nm and/ 
or 248 nm), is normally incident on the AIN:Ag film. Partially, it is absorbed by the film; 
partially, it is reflected and the rest is absorbed by the Si substrate. Figure 5a presents the 
absorption percentage per nm of length for the two wavelengths, where several configura- 
tions of randomly distributed Ag nanoparticles were averaged. The difference in the total film 
absorption, between the 193 and the 248 nm wavelengths, is small, with the 193 nm illumina- 
tion being slightly more absorbed. The different standing wave patterns within the 60 nm film 
are also evident. Of great significance is the distribution of the absorbed light power among 
the different materials (Figure 5b). Considering that the Ag PNPs are dispersed at only 15% 
by volume, it is evident that they are the stronger absorbers in both cases. However, there is 
a significant deviation between the two annealing wavelengths: for 193 nm, a much more bal- 
anced absorption between AIN and Ag occurs, with minimal light energy reaching the Si sub- 
strate. In contrast, for 248 nm, most of the absorption occurs within the Ag inclusions, while at 
the same time, more light energy finds its way to the Si substrate causing significant damage. 
The absorption distribution maps for these two cases are shown in Figure 5c, considering a 
slice of one of the different randomized configurations. The fact that Ag nanoparticles act as 
absorption centres explains why AIN crystallization occurs only within a few nanometres from 
the nanoparticles surface. In addition, the more balanced absorption distribution in the film, as 
well as the minimal absorption by the Si substrate, qualitatively explains the superior annealing 
performance at 193 nm. 
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Figure 5. (a) Absorption distribution for a 60 nm thick AIN:Ag(15%vol)/Si film. (b) Absorption percentages for each 
constituent material within the AIN:Ag/Si structure. (c) Colour coded absorption distribution maps. (d) Absorption 
percentage within each material as a function of nanoparticle size, for 193 nm laser annealing wavelength and (e) Absorption 
percentage within each material as a function of nanoparticle size, for 248 nm laser annealing wavelength. The semi-open 
circles at 20 nm correspond to the nanoparticles having a non-absorbing 2.5 nm crystalline AIN (n = 2.1) layer around them. 


Figures 5d and 5e displays the cumulative absorption in each different material as a function 
of Ag particle size for the two annealing wavelengths. Ordered Ag NP arrays are considered 
here only for simplicity. Evidently, the 193 nm laser exhibits superior performance: unfluc- 
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tuating performance for all particle sizes, more balanced absorption between matrix and 
inclusions, slightly stronger overall film absorption and less substrate absorption. Semi-open 
circles at 20 nm correspond to the nanoparticles having a non-absorbing 2.5 nm crystalline 
AIN (n = 2.1) shell. This results in a small decrease in the matrix absorption due to the higher 
fundamental energy gap of w-AIN. The effect of shell creation thus is towards the direction of 
limiting the film's overall absorption, although it does not seem to be large enough to qualify 
as a self-limiting factor for the annealing process. 


3.1.2.3. The effect of laser fluence 


The data presented in Figure 5 demonstrate that the laser fluence is an essential parameter 
for the LA process. The optical reflectivity spectra of the laser annealed AlN:Ag vs. the laser 
fluence for both processing wavelengths and the corresponding quantitative results are pre- 
sented in Figure 6. Increasing the fluence of the 193 nm wavelength, from 200 to 400 mJ/cm?, 
induces a spectral shift concomitantly with a significant enhancement of plasmon resonance. 
Both effects are respectively associated with changes of the matrix's refractive index and the 
Ag nanoparticles size. The laser fluence decreased the bandwidth (full width at half maxi- 
mum of the fitted LSPR band), indicating a narrower size distribution of Ag nanoparticles due 
to the matrix crystallization, as well as to Ag particle enlargement causing weaker electronic 
confinement. Eventually, at higher fluence (e.g. 600 mJ/cm?), the reflectivity spectrum reveals 
the Si features, indicating that this fluence exceeds the ablation threshold. 


On the other hand, for LA using two pulses at 248 nm, the picture is quite different (Figure 6b). 
At 300 mJ/cm/?, the reflectivity spectrum is similar to the one pulse at 248 nm and 400 mJ/cm? 
(Figure 4c), which is characterized by the absence of LSPR, indicating that the provided energy 
is not enough to cause ripening of Ag and only dissolution is observed. For fluences exceed- 
ing 400 mJ/cm?, a consistent LSPR blue shift is observed due to only a possible reduction of 
the refractive index of the surrounding medium, supporting the Ag outdiffusion to the free 
surface and the formation of Ag nanoparticles surrounded by air. 


3.1.3. Application in optical encoding of information 


The combination of structural modifications and the variation of AIN:Ag plasmonic response, 
introduced by LA, produces strong reflectivity contrast, which has been exploited to encode 
spectrally sensitive patterns [19, 20]. Patterns such as ‘ALN:AG’, ‘ALN’, and ‘AL’ have been 
produced by LA onto films of AIN:Ag, AIN and Al, respectively, deposited either on Si or 
sapphire substrates, using either the 193 nm or the 248 nm beam, and their response (in trans- 
mittance and reflectance) to different light colours was investigated using optical microscopy 
[19, 20]. The reflectivity contrast AR(=R, ,-R ) between treated and untreated areas of an 
AIN:Ag sample demonstrates the level of contrast the human eye can resolve; the spectral 
variations of such AR contrasts for 193 and 248 nm LA are presented in Figure 7b for all cases 
of laser treated AIN:Ag deposited on opaque Si substrates. Depending on the processing 


untreated 


wavelength and LA conditions exceptional flexibility in selecting the required encoding AR 
contrast over the whole visible spectral range can exist. In particular, the spectral selectivity 
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Figure 6. Optical reflectivity spectra with respect to laser fluence for LA at (a) 1 pulse at 193 nm and (b) 2 pulses at 
248 nm. 


is greater when large variations of AR are manifested in the visible spectrum. For example, 
for LA using one pulse of 400 mJ/cm? of 193 nm (Figure 7b, solid blue line), the AR ranges 
from 5% for red (650 nm) to 35% for blue (475 nm) providing very high spectral selectivity. 
These conditions are ideal for spectrally sensitive overt encoding. On the contrary, covert 
encoding requires more subtle spectral contrast AR in order for the pattern to be invisible under 
white light illumination and be revealed only under specific coloured illumination near the 
LSPR region. This can be done by, e.g., LA using two pulses of 500 mJ/cm/? of 248 nm (Figure 7b, 
dashed magenta line). 


The potential of plasmonic encoding has been better demonstrated by optical microscopy 
images, where the sample is illuminated by different coloured light (blue/468 nm, green/532 
nm, red/655 nm) [20]. For this demonstration, AIN:Ag, AIN and Al films were grown on sap- 
phire substrates in order to evaluate both the reflection and transmission reading [20]; the rel- 
evant optical microscopy images are presented in Figure 8. Images A1-D1 demonstrate covert 
encoding in transmission reading, facilitated by LA of an AIN:Ag film at 248 nm. The pattern, 
with a pixel size of 60 um, is almost invisible under white light illumination (A1) and is fully 
uncovered only under blue illumination (D1) where the LSPR is manifested. Images A2-D2 
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Figure 7. (a) Emission spectra of the three red-green-blue (RGB) light sources used for optical microscopy observations. 
(b) Spectral reflectivity contrasts (AR) for laser annealing of AIN:Ag on Si (opaque substrate) at 193 and 248 nm (purple 
and magenta lines, respectively) for various number of pulses and fluences demonstrating the flexibility in selecting the 
required encoding contrast. 


and A3-D3 correspond to overt encoding in both transmission and reflection reading, realized 
by LA at 193 nm. The spectral sensitivity of the plasmonic patterns (A1-D1, A2-D2, A3-D3) 
contrasts the behaviour of patterns encoded by laser ablation in AIN (images A4-D4) and in 
Al (images A5-D5), which do exhibit similar contrasts for all colours, mostly due to scattering 
of light at the pixel edges. In addition to photosensitivity, AIN:Ag(PNPs) films have a unique 
combination of assets such as being harder than steel [20], chemically resistant against organic 
solvents [19, 20], almost having pitchless surface [19, 20] and are resistant to temperatures up 
to 600°C [20, 60]. 


LA@248nm LA@193nm LA@193nm LA@193nm LA@193nm 
In transmission In transmission In reflection . In reflection In reflection 


Figure 8. Optical microscopy images of AIN:Ag, AIN, Alon sapphire using A. white, B. red, C. Green and D. blue illumination, 
in transmission (1, 2) or reflection (3-5) mode on sapphire for 248 nm (1), and 193 nm (2-5) laser annealing, respectively [20]. 
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3.2. Magnetron sputtered AIN/Ag multilayer thin films 


In this section, an alternative pathway is demonstrated, namely, starting from non-plasmonic 
metal/dielectric multilayers, which are then transformed into plasmonic templates by single-shot 
UV LA leading to a subsurface nanoconstruction involving plasmonic Ag nanoparticles embed- 
ded in a hard and inert dielectric matrix. The subsurface encapsulation of Ag nanoparticles pro- 
vides durability and long-term stability. The morphology of the final composite film depends 
primarily on the nanocrystalline character of the dielectric host and its thermal conductivity. The 
emergence of LSPR and its adjustability, depending on the applied fluence, is demonstrated. 


3.2.1. Growth and processing of materials 


Multilayers consisting of alternate thin layers of AIN and Ag were fabricated in an in-house 
built high vacuum system (base pressure 2 x 10 Pa) employing the reactive magnetron sputter- 
ing technique. The AIN crystal structure was varied in the two sets of AIN/Ag by applying two 
different power values to the Al target, 100 and 15 W. This variation lead to different categories 
of sputtered AIN [46]: at 100 W high crystalline wurtzite AIN (w-AIN) was obtained while at 15 
W an amorphous AIN (a-AIN) was obtained. The produced samples were subjected to LA with 
an excimer ArF source. The effect of a single-pulse LA step onto the structural and optical prop- 
erties was investigated as a function of fluence (400-700 mJ/cm’) in an ambient atmosphere. 


3.2.2. Laser processing 


The heating dynamics involved in the multilayer structures during UV LA can be simulated/ 
predicted performing optical and heat transport calculations. UV light propagation was 
solved via the FDTD method [61, 62], and a detailed map of the absorption profile in the 
multilayer and the substrate was obtained. The model system consists of a stratified multi- 
layer structure of alternative dielectric (AIN) and metal layers (Ag) on a Si substrate, capped 
by AIN. The optical response of such a multilayer show high UV absorption (Figure 9a) due 
to silver's inter-band transitions at this spectral range. Under the proper experimental condi- 
tions, this high UV absorption was utilized to locally melt the structure. Due to the immiscibil- 
ity of Ag with AIN, Ag diffusion will drive its transformation into a plasmonic nanocomposite 
with metallic PNPs embedded in a dielectric matrix (Figure 9b). The morphology of the nano- 
composite and its optical response can be tuned by choosing materials and LA parameters 
such as: the nanocrystaline character of the dielectric which is strongly related to its thermal 
conductivity, the metal volume ratio and the total thickness of the multilayer. 


The spatial absorption profile is convoluted with the temporal pulse profile (shown in Figure 11b, 
inset) and solved in the 1D heat transport equation [63]: 


c(z)p(z)0,, T(z, t) = 0,|k(20, T(z, H| + Q(z, t) (1) 


where 0, 0. denote partial derivatives with respect to time (t) and distance (z) from the top 
surface, respectively, and Q(z, t) = a(z) g(t) is the laser heating source term defined by the prod- 
uct of the absorption spatial profile a(z) and the laser pulse temporal profile @(t), with the 
total laser fluence given by: 


f= i odt (2) 
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Figure 9. (a) Spectral absorption of a stratified metal (Ag) dielectric (AIN) structure illustrated in (b) for two metal volume 
ratios 33 and 66%. The high UV absorption can be utilized for the further modification of the nanocomposite. (b) Schematic 
illustration of the laser treatment of the multilayer structure. Under UV laser radiation, the absorbed energy is converted 
into heat, which is diffused by thermal conduction into the Si substrate. The result of the laser treatment is to totally 
transform the multilayer structure into a nanocomposite of embedded plasmonic metallic NPs inside a dielectric matrix. 


c, p and kare the specific heat capacity, the mass density and the thermal conductivity, respectively. 
These three terms have a spatial dependence due to the interchange of different materials. Eq. (1) is 
time-integrated numerically (by a 4th order Runge-Kutta scheme) in order to calculate the explicit 
temperature transient. To do that, we assume that: (i) electronic and lattice temperatures are in 
equilibrium, (ii) radiation and convection losses are negligible, (iii) material properties (absorption, 
heat capacity, mass density and thermal conductivity) remain constant during heating. 


Relevant recent experiments of a single UV laser annealing (LA) on a ceramic/metal (in par- 
ticular AIN/Ag) multilayer film are discussed [17]. This LA step is capable of driving the sub- 
surface modification of the metallic multilayers resulting into LSPR behaviour from PNPs that 
are formed and dispersed in a depth of several nanometre away from the surface. TEM obser- 
vations from the as-deposited samples (Figure 10a and b) verified their multilayer character. 


Figure 10. TEM cross sectional images of: (a) the as-deposited multilayer structure of w-AIN/Ag, (b) the same after LA, 
and (c) the as-deposited multilayer structure of a~-AIN/Ag and (d) the same after LA. LA fluence for all relevant samples 
in this figure was set at 600 mJ/cm?. 
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After LA treatment, the upper part of the multilayer was highly affected. The multilayer struc- 
ture is destroyed and the enlarged NPs are fairly homogenously distributed inside the AIN 
matrix, retaining a high density in the case of a-AIN. On the contrary, in the case of w-AIN 
only, the top of the laser treated part of the sample is affected and is generally depleted of Ag, 
indicating a higher diffusivity of Ag in w-AIN compared to a-AIN. The w-AlN/Ag presents 
a destruction of periodicity affecting only the top six Ag layers, while two more layers were 
partially destroyed. The sample thickness that retained a stratified arrangement is indicated 
by a dashed line in the inset of Figure 10b, which depicts the contrast intensity profile along 
the film thickness. In the case of a-AIN/Ag, this structural re-arrangement is more extensive as 
shown in Figure 10d, and it affected the top 10 Ag layers, while again 2 additional layers were 
partially affected. 


The above considerations make AIN's thermal conductivity a crucial parameter in the 
LA process. k, is sensitive to the deposition method and conditions and can show large 
variations ranging from about 10 W/mK for amorphous AIN [64] to 285 W/mK for single 
crystal [65]. To explore its effect, the temperature transient is plotted (at the top of the mul- 
tilayer) for several different values of the AIN thermal conductivity. A significant difference 
depending on k,n İS found for the peak temperature rise when the laser pulse is on, but it 
quickly diminishes afterwards. The most important effect of k, is shown in the spatial dis- 
tribution of the peak temperature rise plotted in Figure 11(c). In particular, for high values 
of ky the temperature profile is almost constant across the multilayer, even though most 
of the laser fluence is absorbed within the first 100 nm, as shown in Figure 11(a), where 
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Figure 11. (a) UV laser (193 nm) absorption profile in our model structure (21 5nmAg/10nmAIN bilayers on Si), 
(b) normalised temperature transient at the top of the multilayer for various values of k, The inset shows the pulse 
profile of the excimer laser, (c) the normalized peak-temperature spatial distribution for the corresponding cases of 
(b). Large gradients across the film are obtained for low values of k w (d) peak transient temperature distribution 
as a function of distance from the film's top surface, at different incident laser fluence assuming a nonlinear thermal 
conductivity for the multilayer structure (w-AIN and Ag). Solid lines are for the linear case (zero fluence limit) for w-AIN 
(k = 285 W/mK) and for an ultra-low conductance a-AIN (k = 10 W/mK), both for a linear Ag thermal conductivity of 429 
W/mK. Broken lines are for nonlinear w-AIN and Ag for different laser fluence. In the inset we plot the experimental 
temperature dependent thermal conductivities of Ag [50] and w-AlN [48] assumed in this work. 
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the spatial absorption profile within the multilayer is plotted. One can note that for low 
values of kw a significant temperature gradient is developed across the multilayer. Higher 
temperatures close to the surface and lower close to the substrate will result into different 
annealed structures, explaining our experimental results. It should be stressed here that 
up to now only the linear heat transport regime is studied in order to get insight into the 
temperature transients and gradients that can be developed and have not considered the 
thermodynamics of phase changes, e.g. melting and resolidification. Regarding radiation 
and convection losses from the top surface, a simple estimate suffices to show that they are 
of no consequence [20]. 


4, Surface nanostructuring 


Plasmonic templates can be fabricated by LA of a thin metal film on a substrate, resulting 
into nanoparticle arrangements with many plasmonic applications as mentioned in Section 
1 [11, 13, 14, 17, 19, 47-54]. The enabling property of these templates is the coherent surface 
plasmon oscillations excited by the incident light at specific frequencies, with high selectivity 
on the metal nanostructure size and shape. But a selective plasmon excitation is necessarily 
accompanied by selective light absorption, and thus by selective heating of the nanostructure. 
It is natural thus to allow this core property of the target application, to become itself one of 
the drivers in a multi-pulse laser induced self-assembly process. Besides heating, selective 
cooling is also an integral part of restructuring. But as the cooling rate depends on the geo- 
metric features of the nanostructure, a second driver for the self-assembly process is intro- 
duced. The available light absorption mechanisms in the case of nanostructured noble metals 
such as Ag, Au and Cu, can be generally split into (a) inter-band absorption from the metal's 
internal d-electrons, which is typically enabled at UV frequencies and in general is not size- 
dependent, and (b) resonant intra-band absorption, from the conduction electrons due to 
excitation of localized surface plasmon resonances (LSPR), which is strongly size and shape- 
selective. As an example, the overall absorption profile of an Ag nanoparticle (hemispherical 
shape) with diameter of 40 nm on top of a Si substrate is plotted in Figure 12, showing the 
spectral domains of the two absorption mechanisms. Utilization of these absorption modes 
combined with the size selectivity of cooling (~nanostructure surface/volume ratio) consti- 
tutes a toolbox for controlling the self-assembly of nanoparticles. 


If after some laser treatment of a thin silver film on a Si substrate, Ag NPs with specific 
diameters have been formed, what will be the effect of a new incoming laser pulse? For this 
consideration, computer simulations can help explore the potential use of the various modes 
presented above. For computational simplicity a hexagonal lattice of Ag-hemispheres (40 
nm in diameter) on a Si substrate with separation aspect ratios a/d = 1.2 (periodicity/diam- 
eter) have been assumed. First of all, Figure 13a shows the absorption per nm at different 
depths for this structure. The different modes relating to inter-band and intra-band absorp- 
tions in the nanoparticle, as well as the absorption in the substrate, are clearly shown. The 
absorption spatial distribution is then time-integrated with the appropriate laser temporal 
profile for each wavelength, depicted in Figure 13c. The experimentally acquired laser pulse 
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Figure 12. Simulated absorption spectrum of a 40 nm diameter Ag nanoparticle on top of a Si substrate. Inter-band 
absorption dominates the red shaded region, while LSPR driven absorption is dominating the green shaded regions (at 
~350 nm the quadrupole 428 and ~580 nm the dipole contributions, respectively). 


temporal profile from an excimer source was used for the 150-180nm range whereas the 
pulse temporal profile of a Nd:YAG source was used for the 300-900 nm range in a numeri- 
cal solution of the numerically 1D heat diffusion equation presented in Eq. (1). Thus, we 
get the temperature transient profile at every point in the structure in degrees/mJ of laser 
fluence. In Figure 13d, the temperature transients for the two primary wavelengths used in 
recent experimental studies are presented. This way the mapping of the maximum tempera- 
ture rise as a function of particle size, laser wavelength and laser fluence can be constructed. 
Bear in mind that the temperature transient follows the actual pulse temporal shape, but a 
significant transient tail also persists for longer times. The above discussion indicates that 
the processes discussed here are power and fluence dependent respectively, making the 
explicit time-dependent simulations essential. For a given laser fluence, direct information 
regarding particle heating can be extracted, and the probability of a nanoparticle melting or 
not can be assessed. 


The global picture can then be revealed by considering how nanoparticle arrays of differ- 
ent sizes of nanoparticles behave to different wavelength illuminations. This is plotted in 
Figure 13e, showing the peak transient temperature, per unit of laser fluence, developed in 
the nanoparticle as a function of both the wavelength and the nanoparticle size. The dipole 
contribution is weaker than the contributions at shorter wavelengths (quadrupole LSPR and 
Si inter-band absorption), however, within a specific annealing plan one can always irradiate 
at one single wavelength, thus the strength of each band is of no importance other than defin- 
ing the fluence that each laser wavelength should be used. Upon melting and resolidification, 
a random distribution of sizes will emerge. However, repeated laser pulses at a particular 
wavelength will lead to the melting of one particular size group and thus to its eventual 
extinction. The main concept behind this approach is, thus, repeated laser treatment, target- 
ing different size groups with different laser wavelengths, driving the final size distribution 
towards a predetermined one. This process has been termed as ‘modification of nanopar- 
ticles arrays by laser-induced self-assembly’ (MONA-LISA) [47] and has been experimentally 
demonstrated. 
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Figure 13. (a) Vertical absorption (% per nm) profile for the 40 nm diameter Ag nanoparticle. Indicated with the two 
dashed lines are the two primary wavelengths used (248 and 532 nm). (b) Absorption percentage per nm profile in the 
model structure for the two wavelengths used. (c) Temporal profile of the two different laser sources. (d) Corresponding 
temperature transients developed in the Ag NP (normalized to the laser fluence). (e) Global map showing the peak 
temperature developed in an Ag NP on top of Si versus wavelength and NP size. 


As a first step, a single shot UV LA process was utilised, to initialize nanoparticle formation 
from a semi-continuous metal thin film, by laser induced self-assembly (LISA). This resulted 
into the restructuring of the Ag thin film into nanoparticle arrangements in macroscopic scales 
(several mm’). The LISA process with a UV laser (248 nm) is strongly fluence-dependent as 
demonstrated by the SEM images of Figure 14a. The Ag NP size distributions are generally 
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Figure 14. (a) SEM images of samples processed with LISA (1 pulse, 248 nm LA) and UV MONA-LISA (10 pulses, 248 
nm) under various applied fluences. (b) The effect of LISA and UV MONA-LISA on the optical reflectivity spectra 
of a 10 nm Ag thin film under various applied fluences. (c) Surface coverage of the nanoparticles in relation to the 
particle diameter for all the cases of (a). (d) Critical nanoparticle diameter (the diameter value under which every 
particle remains practically cold during the UV LA process, hence these particles cannot be further manipulated 
upon irradiation) as a function of applied laser fluence. Inset depicts a portrait of Da Vinci's Mona Lisa: a microscopic 
pattern engraved on silicon wafer in a single pixel processing using series of photolithographic masks. (e) A digital 
photo of a grid of LISA / UV MONA-LISA spots of varying fluence (x-axis) and number of pulses (y-axis) produced 
ona 10 nm Ag thin film. 


bimodal, consisting of one set of coarse and one set of fine NPs. Increasing the laser fluence from 
350 to 800 mJ/cm? resulted into a decreasing difference between the two NP sets, as shown by 
the quantified SEM images (Figure 14c, purple, cyan, yellow and red lines). The modification 
of NP assemblies by multiple laser pulses (MONA-LISA process) involves the utilization of the 
selective optical absorption and the heat dissipation, as presented in the theoretical approach 
above. For this consideration, the initial assemblies of NPs were subjected to a series of UV and/ 
or VIS LA steps. 


The case of UV MONA-LISA is first considered, i.e. processing by multiple UV (248 nm) 
laser pulses, and the effect of inter-band absorption is correlated with the size-selective 
heat dissipation. Figure 14b shows the optical reflectivity spectra of a 10 nm Ag film 
recorded after UV MONA-LISA applying 10 successive UV laser pulses at various fluences 
(magenta, blue, green and orange lines for 350, 500, 650 and 800 mJ/cm’, respectively). 
Comparing these spectra with those from the previous LISA processing revealed that 
the consecutive pulses transform the two reflectivity peaks into one, indicating a gradual 
transformation of the bimodal distributions into unimodal ones. This was also confirmed 
by the corresponding SEM images and the size distribution histograms in Figure 14a and c 
(magenta, blue, green and orange bars), respectively. The reshaping of the particles is pre- 
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dominantly due to the disappearance of the larger particles. The physics behind that obser- 
vation is based on the stronger heating of the larger NPs and faster heat dissipation from 
the finer NPs due to their higher effective surface (surface to volume ratio). As a result, the 
finer nanoparticles remain cold during UV MONA-LISA. To quantify this, photothermal 
calculations [47] were employed, in order to determine the equitemperature contour lines 
vs. nanoparticle diameter, from which the critical size under which the particles remain 
practically cold during the UV MONA-LISA. Considering the bulk Ag melting point of 
960°C, the calculated values of the critical diameter vs. the laser fluence (248 nm) are dis- 
played in Figure 14d. We remark that the critical diameter is highly dependent on the 
thermal conductivity of the substrate; less thermally conductive substrate would hinder 
heat dissipation. The UV MONA-LISA process is thus based on the ability of the larger par- 
ticles to raise their temperature above their melting point, and consequently re-solidifying 
forming smaller particles with a stochastic size distribution, while the smaller ones remain 
below the melting threshold and therefore remain practically unaffected. As a result, suc- 
cessive pulses suppress the population of the coarse particles, gradually transforming 
them to finer particles via stochastic size redistribution. The wide range of size distribu- 
tions achieved by LISA and UV MONA-LISA results in a vast colour palette of plasmonic 
responses as those presented in Figure 14e (real visual appearance of the Ag nanoparticle 
assemblies). Combining LISA and MONA-LISA process with projection masks enables the 
production of multicolour images of high spatial resolution and colour contrast. Such an 
example is presented in the inset of Figure 14d, depicting the portrait of Da-Vinci's famous 
‘Mona Lisa’. 


The implementation of MONA-LISA employing exclusively UV pulses has the inherent lim- 
itation of the UV absorption being rather insensitive to nanoparticle size, and as a result, the 
size selectivity of the process is based exclusively on the size selectivity of the heat dissipa- 
tion. Thus, only the large particles are affected by the successive laser pulses. Finer particles 
can be further manipulated by exploiting their selective optical absorption via LSPR using 
a visible wavelength, resulting in larger optical absorption and hence higher temperature 
rise, as predicted by the computational results presented in Figure 13. The optimal condition 
of narrowing the LSPR bandwidth at 532 nm was achieved by applying 500 pulses at 125 
mjJ/ cm’, which shows that the VIS MONA-LISA is slower compared to the UV MONA-LISA 
that reached at steady state at 10 pulses. This is well expected considering the computational 
results of Figure 13, where a lower temperature rise is evident at visible wavelengths com- 
pared to UV processing. 


Exploring the VIS MONA-LISA one can start from a UV MONA-LISA process (248 nm, 2 
pulses at 350 mJ/cm?) which demonstrates a double reflectance peak (due to bimodal nano- 
structuring) notably with no particular response close to 532 nm, as shown in Figure 15a 
(black line). Applying, VIS MONA-LISA to this area, with 500 pulses (i), 2000 pulses (ii) and 
4000 pulses (iii) results in gradually re-adjusting the bimodal distribution to one with a peak 
very close to the laser processing wavelength (532 nm). The corresponding SEM images are 
shown in Ref. [47]. Analysis of these SEM images results in the histograms of Figure 15b. It is 
noteworthy that the VIS MONA-LISA probes particularly the finer Ag nanoparticles (<20 nm) 
confirming the computational analysis presented in Figure 13. 
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Figure 15. (a) Optical reflectivity spectra of LISA (2 pulse, 248 nm at 350 mJ/cm? - black line), used as a template for the 
subsequent Visible MONA-LISA (500 pulses (red line), 2000 pulses (blue line) and 4000 pulses (magenta line), 532 nm 
at 125 mJ/cm’). The green dashed line sets the 532 nm mark. (b) Surface coverage of the nanoparticles in relation to the 
particle diameter for each laser treatment of (a). 


5. Conclusions 


Functional plasmonic templates can be fabricated by LA. In particular, two distinct cases have 
been discussed: laser nanostructuring of either ceramic/metal multilayers and composites or 
ultra-thin metal films results on sub-surface or surface plasmonic patterns, respectively, with 
many potential applications. 


Theoretical investigations have highlighted the underline physics of laser-induced sub-surface 
and surface plasmonic nanostructuring: the absorption profile of the film and the temperature 
spatial distribution strongly depend on the structure parameters such as the metal volume 
ratio, the total thickness of the film and the thermal conductivity of the dielectric. By properly 
designing these parameters, combined with tuning the LA conditions (wavelength, number 
of pulses and fluence) one can arrive at plasmonic templates with predetermined morphol- 
ogy and optical response. The latter was demonstrated in experiments in which either AIN:Ag 
nanocomposites or multilayers of alternative AIN/Ag layers were subjected in UV LA leading 
to a sub-surface formation of PNPs sealed in an AIN matrix with a localized surface plasmon 
resonance response. 


An extensive experimental and theoretical investigation has highlighted the photothermal 
processes involved in laser-induced surface plasmonic nanostructuring. Specifically, tem- 
plates consisting of hemi-spherical Ag PNPs on Si substrate, with different arrangements 
and sizes, for both UV and optical irradiation have been explored. Different absorption 
mechanisms, which are size selective and enable geometric specificity, have been identified: 
(a) inter-band absorption taking place at the UV zone, probing the metal's internal d-elec- 
trons and (b) intra-band absorption taking place at the near-UV, and visible zone, probing 
the metal's free electrons at the specific spectral position of LSPR. In particular, by sequen- 
tially tuning the laser wavelength into resonance with different physical absorptions, we 
can selectively target the melting and re-solidification of different particle size groups. This 


process was termed: “modification of nanoparticle arrays by laser-induced self-assembly” 
(MONA-LISA). 
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